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Abstract — In this paper, we study the performance of blind 
adaptive beamforming algorithms for smart antennas in realistic 
environments. A constrained constant modulus (CCM) design 
criterion is described and used for deriving a recursive least 
squares (RLS) type optimization algorithm. Furthermore, two 
kinds of scenarios are considered in the paper for analyz- 
ing its performance. Simulations are performed to compare 
the performance of the proposed method to other well-known 
methods for blind adaptive beamforming. Results indicate that 
the proposed method has a significant faster convergence rate, 
better robustness to changeable environments and better tracking 
capability. 
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I. Introduction 

Recently, adaptive beamforming has attracted considerable 
attentions and been widely used in wireless communications, 
radar, sonar, medical imaging and other areas [Q]. Many 
existing methods have been presented in different communi- 
cation systems (2)-|[5]]. Blind adaptive beamforming, which is 
intended to form the array direction response without knowing 
users' information beforehand, is an encouraging topic that 
deals with interference cancellation, tracking improvement and 
complexity reduction. 

The linearly constrained minimum variance (LCMV) 
method, with multiple linear constraints [6], is a common 
approach to minimize the beamformer output power while 
keeping the signal of interest (SOI) from a given direction. 
However, because of the required input data covariance matrix, 
the LCMV beamformer cannot avoid complicated compu- 
tations, especially for large input data and/or large sensor 
elements. Also, this method suffers from slow convergence 
due to the correlated nature of the input signal [?]. 

Choi and Shim [9| proposed another computationally effi- 
cient algorithm based on the stochastic gradient (SG) method 
for finding the optimal weight vector and avoiding the esti- 
mation of the input covariance matrix. As shown in ||9l, a 
cost function is optimized according to the minimum variance 
subject to a constraint that avoids the cancellation of the SOI, 
i.e. the so-called constrained minimum variance (CMV). Nev- 
ertheless, this algorithm still cannot avoid slow convergence 
rate. Furthermore, another noticeable problem is how to define 
the range of step size values. The small value of step size will 



lead to slow convergence rate, whereas a large one will lead 
to high misadjustment or even instability. 

In 1 10 1 SG algorithms using the constant modulus (CM) cost 
function are reviewed by Johnson et al. for blind parameter 
estimation in equalization applications. Similarly, the CM 
approach exploits the low modulus fluctuation exhibited by 
communications signals using constant modulus constellations 
to extract them from the array output. Although it adapts 
the array weights efficiently regardless of the array geometry, 
knowledge of the array manifold or the noise covariance 
matrix, the CM-based beamformer is quite sensitive to the step 
size. In addition, the CM cost function may have local minima, 
and so doesn't have closed-form solutions. Xu and Liu ifTTI 
developed a SG algorithm based on constrained constant 
modulus (CCM) technique to sort out the local minimum 
problem and obtain the global minima. But they still cannot 
find a satisfied solution in terms of the slow convergence. 

To accelerate convergence, Xu and Tsatsanis [12 | employed 
CMV with the recursive least squares (RLS) optimization 
technique and derived the CMV-RLS algorithm. It turns out 
that this method exhibits improved performance and enjoys 
fast convergence rate. However, it is known to experience per- 
formance degradation if some of the underlying assumptions 
are not verified due to environmental effects. The signature 
mismatch phenomenon is one of these problems. 

In this work, we propose a constrained constant modulus re- 
cursive least squares (CCM-RLS) algorithm for blind adaptive 
beamforming. The scheme optimizes a cost function based on 
the constant modulus criterion for the array weight adaptation. 
We then derive an RLS-Type algorithm that possesses better 
performance than those of previous methods. We carry out a 
comparative analysis of existing techniques and consider two 
practical scenarios for assessment. Specifically, we compare 
the proposed method with the existing CMV-SG (9), CMV- 
RLS Oil and CCM-SG fl4l . 

The remaining of this paper is organized as follows. In the 
next section, we present a system model for smart antennas. 
Based on this model, we describe the existing SG algorithms 
based on the constrained optimization of the minimum vari- 
ance and constant modulus cost functions in Section III. In 
Section IV, the RLS-type algorithms, including the proposed 
algorithm, are derived. Simulation results are provided in 



Section V, and some conclusions are drawn in Section VI. 
II. SYSTEM MODEL 

In order to describe the system model, let us make two sim- 
plifying assumptions for the transmitter and receiver models 
fl6l . First, the propagating signals are assumed to be produced 
by point sources; that is, the size of the source is small with 
respect to the distance between the source and the sensors 
that measure the signal. Second, the sources are assumed to 
be in the "far field," namely, at a large distance from the 
sensor array, so that the spherically propagating wave can 
be reasonably approximated with a plane wave. Besides, we 
assume a lossless, nondispersive propagation medium, i.e., a 
medium that does not attenuate the propagating signal further 
and the propagation speed is uniform so that the wave travels 
smoothly. 
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Fig. 1. Adaptive beamforming structure for ULA. 

Let us consider the adaptive beamforming scheme in Fig. 
Q] and suppose that q narrowband signals impinge on the 
uniform linear array (ULA) of rn (q < m) sensor elements 
from the sources with unknown directions of arrival (DOAs) 
6*o,- . . ,8q-i- The ith snapshot's vector of sensor array outputs 
can be modeled as ITTll 

x(i) = A(6)s(i) + n(i), i=l,...,N (1) 

where 6 = [9 , . . . , 9 q -i] T G C qxl is the vector of the 
unknown signal DOAs, A(9) = [a{6 ), a(^_i)] e C mxq 
is the complex matrix composed of the signal direction 
vectors a((9 4 ) = [i^-^i"*,.-^" 2 "^ 1 ^™ 9 '] 5, G 
C mxl , (i = 0, . . . , q — 1), where A is the wavelength and d = 
A/2 is the inter-element distance of the ULA, s(i) G lZ qxl is 
the real value vector of the source data, n(i) G C mxl is the 
complex vector of white sensor noise, which is assumed to be 
a zero-mean spatially and white Gaussian process, N is the 
number of snapshots, and (-) T stands for the transpose. The 
output of a narrowband beamformer is given by 

y(t) = w H x{i) (2) 



where w — [wi, . . . , w m ] T G C mxl is the complex weight 
vector, and (-) H stands for the Hermitian transpose. 

Constrained adaptive beamformers are the ones that min- 
imizes the objective function J w subject to a set of linear 
constraints. That is 

w = arg min J w subject to C H w = f (3) 

w 

where w as mentioned before is the vector of coefficients, 
C G C mxq is a constraint matrix, and / G C qxl is a constraint 
vector with respect to sources. 

III. Stochastic Gradient Algorithms 

The purpose of SG algorithms is to get an acceptable output 
performance and reduce the complexity load by avoiding the 
correlation matrix estimation and inversion. We first describe 
the CMV and CCM algorithms on the basis of the SG method, 
which are called CMV-SG and CCM-SG, respectively. These 
two methods employ different cost functions, which are pre- 
sented below. 

A. CMV-SG algorithm 

The CMV-SG algorithm chooses the beamformer to mini- 
mize the sum of the weighted array output power. The con- 
strained optimization means that the technique minimizes the 
contribution of the undesired interferences while maintaining 
the gain along the look direction to be constant. For elaborating 
these principles clearly, we set C = a(9o), where 6q is the 
DOA of the desired signal, and then the received / reduces to 
a constant value, say 1. The expected result of the algorithm is 
to update the weight vector with respect to the input data for 
minimizing the cost function of the algorithm. So, in the SG 
approach, the optimization problem of the CMV is represented 
as follows according to the cost function 

J = w H (i)x(i)x H (i)w(i), i = l,...,N ^ 

subject to w H (i)a(Oo) = 1 

Using the method of Lagrange multipliers, the solution can 
be obtained by setting the gradient term of (4) with respect to 
w(i) equals to zero and using the constraint w H (i)a(9o) = 1 

w(i + 1) = w(i) - w*(i)[x(i) - a H (d Q )x(i)a(9 )} (5) 

where /i is the step size value of the CMV-SG algorithm and 
* denotes complex conjugate. 

B. CCM-SG algorithm 

For the CCM-SG algorithm, we consider the cost function 
as the expected deviation of the squared modulus of the array 
output to a constant, say 1 . Compared with the CMV method, 
the CCM cost function is simply a positive measure of the 
average amount that the beamformer output deviates from 
the unit modulus condition Q. By using the same constraint 
condition, the cost function of CCM-SG can be expressed as 

J = (|^)| 2 -l) 2 , i = l,...,N 
subject to w H (i)a(9o) = 1 



Also, using the same operation as that of CMV-SG, the 
update equation of weight vector can be obtained 

w(i + l)=w(i)-»(\y(i)\ 2 -l) y *(i) 
■[x(i)-a H (9 Q )x(i)a(9 Q )] 

It will be shown via simulations that, by employing the 
constant modulus property and the constrained condition, the 
performance of CCM-SG algorithm is better than that of CMV- 
SG one. 

The SG algorithms employ stochastic inputs instead of 
deterministic inputs in the steepest descent method and do 
not use the desired response as the least mean square (LMS) 
method. It is well known that the algorithms based on the SG 
technique generally have low computational complexity. How- 
ever, both the slow convergence behavior and the sensitivity 
to the step size are drawbacks. 

IV. Recursive Least Squares (RLS) Type 
Algorithms 

To overcome the shortcomings of the SG algorithms, we 
propose a fast converging RLS-type algorithm using the CCM 
criterion for the array weight adaptation. An important feature 
of this method is that its rate of convergence is typically 
an order of magnitude faster than those of the simple SG 
algorithms, due to the fact that RLS technique whitens the 
input data by using the inverse correlation matrix of the data. 
Despite the initial presentation of the weight update includes 
the correlation matrix inversion, recursive computations can be 
employed to update the value of the correlation matrix of the 
input data step by step. Then, according to the matrix inversion 
lemma defined in |18|, we obtain the recursive equation for 
the inverse of the correlation matrix and this reduces by one 
order of magnitude the high computational burden introduced 
by matrix inversion. 

A. CMV-RLS algorithm 

Here, we extend the least squares (LS) method to develop 
RLS-Type algorithm for solving the ill-posed problem [18|. 
In this section, we describe the CMV-RLS algorithm first 
and then derive the proposed CCM-RLS method by using the 
same idea. The LS-type cost function of CMV algorithm is 
described by 

i 

J = Y^a i - l w B {i)x{l)x H {l)w(i), i = l,...,N 

1=1 W 
subject to w H (i)a(9a) = 1 

where a is a positive parameter close to, but less than, unity. 

Equation (8) will lead to the LS algorithm. Using the 
method of Lagrange multipliers, setting the gradient term of 
cost function with respect to w(i), making it equal to zero and 
using the constraint w H (i)a(9o) = 1, we obtain the solution 

w(i) = [a H (9 )R- 1 (i)a(9 Q )]- 1 R- 1 (i)a(9 Q ) (9) 
where R(i) = £} =1 a?' 1 x{l)x H {I) 6 C mxm . 



Following the recursion for updating the correlation matrix, 
it can be performed as 

R(i) = aR{i - 1) + x(i)x H {i) (10) 

Then, as mentioned before, can be estimated by 

exploiting a relation in matrix algebra as the matrix inversion 
lemma. If we define P(i) = the RLS solution of 

CMV-RLS is available. 

w(i) = [a H (9o)P(t)a(9o)}- 1 P(t)a(9 ) (11) 

Both Xu and Tsatsanis have proved that the CMV-RLS 
method, compared with SG algorithms, exhibits better effects 
on both performance and convergence behavior in 11121 . How- 
ever, these properties are adversely affected by dynamic envi- 
ronments in practice, which further complicated the estimate 
of the statistics. 

B. Proposed CCM-RLS algorithm 

In order to conquer the problem occurred in CMV-RLS, 
we propose an improved algorithm which uses the constant 
modulus property. 

Consider the optimization equation shown below 

i 

J = ^ a i -'[K(i) a; (0| 2 -l] 2 , i = l,...,N 
i=i 

subject to w H (i)a(8o) = 1 

where a is similar to that in the CMV-RLS algorithm. Here, 
we do not need to bother with the sensitivity of step size value 
since the cost function above even does not contain this term. 
Also, it is easy to choose a due to its small range. As can be 
seen from the CCM-RLS cost function, a closed-form solution 
(global minimum) is not possible because it is a fourth-order 
function with a more complicated structure, i.e. multiple local 
minima. The constraint condition used in (12) ensures that we 
can get the optimal solution. 

Equation (12) is related to the LS algorithms. Using the 
same operation as that of CMV-RLS, we can get the conven- 
tional solution of our proposed method 

w(i) = [a H (9 )R- 1 (t)a(9oT 1 R- 1 (t)a(9 a ) (13) 

where R(i) = 2J2] =1 a l ^ l [\w H (i)x{l)\ 2 - l]x(l)x H {I) e 
C mxm . LS estimation, like the method of least squares, is an 
ill-posed inverse problem. To make it "well posed," we need 
to renew the correlation matrix R(i) [18] as follows 

i 

R(i) = 2^a i - l [\w H {i)x{l)\ 2 -l]x(l)x H (l) + 6a i I. (14) 

2 = 1 

where 6 is a positive real number called the regularization 
parameter and J is the mxm identity matrix. The second term 
on the right side of (14) is included to stabilize the solution to 
the algorithm by smoothing the solution and has the effect of 



TABLE I 
CCM-RLS ALGORITHM 



Initialization 


w(0) = 

P(0) = 8- 1 I mxm 

_ / small positive constant for high SNR 
— \ large positive constant for low SNR 


Update 

(For each 
instant of time 
i = 1, . . . , N) 


7r(i) = P(i - l)x(i) 
e(i) = [\w H (i)x(i)\ 2 - 1] 

*W ~ a/2e(i) + x H (i)ir{i) 

P(i) = a ~ 1 P(i - 1) - a- 1 k(i)x H {i)P(i - 1) 
w(l) = [a ff (e )P(i)a(e )]- 1 P(i)a(e ) 



making the correlation matrix R{i) nonsingular at all stages 
of the computation. 

Let e(i) — [\w H (i)x(l)\ 2 — 1] and follow the recursion for 
updating the correlation matrix, (14) can be expressed as 

R(i) = aR(i - 1) + 2e(i)x(i)x H (i) (15) 

By using the matrix inversion lemma in (15), we can obtain 
the inverse of R(i) 

R- 1 (i) = a- 1 R- 1 (i-l) 

_ a-iR^ji - l)z(i)z fl (t)JT 1 (i - 1) (16) 
l/2e(i) + a- 1 x H {i)R~ 1 {i - l)x(i) 

Here, we still define P(i) = For convenience of 

computation, we also introduce a vector k(i) £ C mxl as 



fe(i) 



a _1 P(z - l)x(i) 



l/2e(i) + a- 1 x H (i)P(i - l)x(i) 
Using these definitions, we may write (16) as 

P(i) = a~ x P(i - 1) - a^ktyx 11 (i)P(i - 1) 



(17) 



(18) 



Until now, we develop a recursive equation to update the 
matrix P(i) by incrementing its old value. Finally, using the 
fact that P(i) equals to we get the proposed RLS 

solution 



w(i) = [a H (9 )P(i)a(0 Q )}- 1 P( i )a(9 Q ) 



(19) 



Equations (17)-(19), collectively and in that order, constitute 
the derived CCM-RLS algorithm, as summarized in Table U 

We note that, (17) and (18) enable us to update the value of 
the vector k(i) itself. An important feature of this algorithm 
described by these equations is that the inversion of the 
correlation matrix R(i) is replaced at each step by a simple 
scalar division. Also, in the summary presented in Table [I] the 
calculation of the vector k(i) proceeds in two stages: 

* First, an intermediate quantity, denoted by 7r(i), is com- 
puted. 

> Second, is used to compute k(i). 

This two-stage computation of k(i) is preferred over the 
direct computation of it using (16) from a finite-precision 
arithmetic point of view 11181 . 



To initialize the CCM-RLS method, we need to specify two 
quantities: 

• The initial weight vector w(Q). The customary practice 

is to set w(0) = 0. 
« The initial correlation matrix R(0), Setting i — in (14), 

we find that R(0) = SI. 
In terms of complexity, the LS requires (D(m 3 ) arithmetic 
operations, whereas the RLS requires 0(m 2 ). Furthermore, 
we can notice that the step size fi in the SG algorithms is 
replaced by R 1 (i). This modification has a significant impact 
on improving the convergence behavior. 

V. Simulations 

The performance of the proposed CCM-RLS algorithm 
is compared with three existing algorithms, namely CMV- 
SG, CCM-SG and CMV-RLS, in terms of output signal-to- 
interference-plus-noise ratio (SINR), which is defined as 



SINR(i) 



w H (i)R s {i)w{i) 
w H (i)R i+n (i)w(i) 



(20) 



where R s (i) is the autocorrelation matrix of the desired signal 
and Ri +n (i) is the cross-correlation matrix of the interference 
and noise in the environment. 

An ULA containing m = 16 sensor elements with half- 
wavelength spacing is considered. The noise is spatially and 
temporally white Gaussian noise with power = 0.01. 
For each scenario, K = 100 iterations are used to get each 
simulated point. In all simulations, the DOA of the SOI is 
6>o = 20° and the desired signal power is ctq = 1. The 
interference-to-noise ratio (INR), in all examples, is equal to 
lOdB. The BPSK scheme is employed to modulate the signals. 
The value of a was set equal to 0.998 in order to optimize 
the performance of the RLS-type algorithms. 

Fig. [2] includes two experiments. Fig. |2ja) shows the output 
SINR of each method versus the number of snapshots, whose 
total is 1000 samples. There are two interferers with DOAs 
0\ = 40° and 82 = 60°. In this environment, the actual 
spatial signature of the signal is known exactly. The result 
shows that the RLS-type algorithms converge faster and have 
better performances than those of the SG algorithms. The 
steering vector mismatch scenario is shown in Fig. |2jb). We 
assume that this steering vector error problem is caused by 
look direction mismatch The assumed DOA of the SOI 
is a random value located around the actual direction, whose 
mismatch is limited in a range of 1°. Compared with Fig.|2ja), 
Fig-Eb) indicates that the mismatch problem leads to a worse 
performance for all the solutions. The CMV-RLS method 
is more sensitive to this environment, whereas the proposed 
CCM-RLS algorithm is more robust to this mismatch. 

In Fig. [3] the scenario is the same as that in Fig. [2 a) for the 
first 1000 samples. Two more users, whose DOAs are 83 = 30° 
and 84 = 50°, enter the system in the second 1000 samples. As 
can be seen from the figure, SINRs of both the SG and RLS- 
Type algorithms reduce at the same time. It is clear that the 
performance degradation of the proposed CCM-RLS is much 
less significant than those of the other methods. In addition, 



RLS-type methods can quickly track the change and recover 
to a steady-state. At 2000 samples, two interferers with DOAs 
65 = 25° and 6*6 = 35° enter the system whereas one interferer 
with DOA 64 = 50° leaves it. The simulation shows nearly 
the same performance as that of the second stage. It is evident 
that the output SINR of our proposed algorithm is superior to 
the existing techniques. This figure illustrates that the CCM- 
RLS is still better after an abrupt change, in a non-stationary 
environment where the number of users/interferers suddenly 
changes in the system. 



(a) SNR=20dB, q=3 users, m=16 



(b) SNR=20dB, q=3 users, m=16, mismatch=1 
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Fig. 2. Output SINR versus number of snapshots for (a) ideal steering vector 
condition (b) steering vector with mismatch. 



SNR=20dB,q =3 users(1-1000); 
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Fig. 3. Output SINR in a scenario where additional interferers suddenly 
enter and/or leave the system. 



VI. Conclusions 

In this paper, a new algorithm enabling blind adaptive beam- 
forming has been presented to enhance the performance and 



improve the convergence property of the previously proposed 
adaptive methods. Following the CCM criterion, a RLS-type 
optimization algorithm is derived. In the place of step size, we 
employ the correlation matrix inversion instead for increasing 
the convergence rate. Then, matrix inversion lemma was used 
to solve this inversion problem with reduced complexity. We 
considered different scenarios to compare CCM-RLS algo- 
rithm with several existing algorithms. Comparative simulation 
experiments were conducted to investigate the output SINR. 
The performance of our new method was shown to be superior 
to those of others, both in terms of convergence rate and 
performance under sudden change in the signal environment. 
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